Dynamics of the (spin-) Hall effect in topological insulators and graphene 
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A single two-dimensional Dirac cone with a mass gap produces a quantized (spin-) Hall step in the 
absence of magnetic field. What happens in strong electric fields? This question is investigated by 
analyzing time evolution and dynamics of the (spin-) Hall effect. After switching on a longitudinal 
electric field, a stationary Hall current is reached through damped oscillations. The Hall conductivity 
remains quantized as long as the electric field (E) is too weak to induce Landau-Zener transitions, 
but quantization breaks down for strong fields and the conductivity decreases as l/\f~E. These 
apply to the (spin-) Hall conductivity of graphene and the Hall and magnetoelectric response of 
topological insulators. 

PACS numbers: 73.20.-r,72.80.Vp,73.43.-f 



The unique electronic properties of graphene can be 
traced back to the the pseudo-relativistic Dirac equation 
and its linear energy dispersion with zero bandgap. It ex- 
hibits a plethora of interesting and fascinating physical 
phenomena related to electric and heat transport, mag- 
netic field effects, valley and spintronics[l[. The 'half- 
integer' quantum Hall effect, in spite of its environmental 
fragility, has been observed at room temperature due 
to the unusual Landau quantization of Dirac electrons. 
With spin-orbit coupling taken into account, graphene 
in principle realizes a spin-Hall insulator belonging to 
the class of topological insulators (TI). 

Moreover, Dirac electrons also occur as surface states 
of three-dimensional TI0-Q, close relatives of integer 
quantum Hall states. These materials are predicted to 
display a variety of peculiar phenomena, such as spin- 
and surface quantum Hall effects and the closely re- 
lated topological magnetoelectric effect Q, allowing for 
the control of magnetization by electric field. As op- 
posed to the even number of Dirac cones in graphene, 
three dimensional TI can have an odd number of Dirac 
cones on a surface. Due to time reversal symmetry, these 
states are robust with respect to non-magnetic disorder, 
similarly to how pair breaking in s-wave superconductors 
is prohibited by potential scatterers. 

The hallmark of (pseudo-) relativistic massive Dirac 
electrons Q is a single quantum Hall step around half- 
filling in the absence of magnetic field between o~ xy = 
±e 2 /2h. We ask how this picture gets modified in the 
presence of strong electric fields. Generally, a driving 
electric field can produce a sizable density of electron 
and hole excitations around the Dirac point in a highly 
non-thermal, non-stationary momentum distribution[9(. 
Consequently, the longitudinal transport of Dirac elec- 
trons features Klein-tunneling [To| and Schwinger's pair 
production^, in a stationary or time-dependent 

framework, when the electric field is represented by 
a static scalar or a time-dependent vector potential, 
respectively. The latter approach directly yields the 



non-equilibrium momentum distribution and the time- 
dependent current at finite electric fields. While it does 
not use any kind of equilibrium or out of equilibrium re- 
sponse formalism (Kubo/Landauer), it still reproduces 
known results and makes predictions for the non-linear 
behaviour of the electric current as an example P.fl4l|. On 
the other hand, a strong electric field alters not only the 
longitudinal transport [9(, but is expected to modify the 
transverse conductivity, involving the non-equilibrium 
quantum (spin-) Hall breakdown. Common wisdom tells 
that while there are no power law corrections to the in- 
teger Hall conductivity for weak electric fields, with its 
quantization 'topologically' protected, there can be expo- 
nentially small corrections. When these grow with field, 
quantization breaks down. 

To consider the problem in detail, we elaborate on 
the time evolution of the Hall current for massive Dirac 
Fcrmions, after switching on a longitudinal electric field. 
We show that a stationary transverse current develops 
for long times, characterized by a quantized Hall con- 
ductivity for weak fields, crossing over to a strongly field 
dependent Hall response with increasing field. This re- 
sult applies to the quantum (spin-) Hall breakdown of 
graphene [3] as well as for the related Q surface Hall and 
magnetoelectric effect in TI. 

The low energy description around the K point in the 
Brillouin zone of graphene Q or on the surface state of a 
3D TI0, El (after a tt/2 rotation of the spin around z), 
in the presence of a uniform, constant electric field (E > 
0) in the x direction [switched on at t = 0, through a 
time dependent vector potential A(t) = (A(t), 0,0) with 
A(t) = EtQ(t)] is written as 

ihd t ^p(t) = (v F (p x - eA(t)), v FPyi A) • tr%(t), (1) 

where vf is the Fermi velocity, and the Pauli matrices 
(cr) encode the two sublattices[l[ of the honeycomb lat- 
tice in graphene, or the physical spin in TI. A > is 
the mass gap, originating from the intrinsic spin-orbit 
coupling (SOC) in grapheneQ, or from a thin ferromag- 
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netic film covering the surface of TI, lifting the Kramer's 
degeneracy of the Dirac point. 

To make our analysis more transparent, we per- 
form a two-step unitary transformation, U = XJ\Ui- 
Firstly, a time independent rotation around the a x 
axis as U\ = C+ — i<r x C-, with C± = [1 ± 
[vfVvI \J ( v FP y ) 2 + A 2 )] 1 / 2 /^/2, and secondly a time 
dependent one, bringing us to the adiabatic basis 
as U2 — exp(—itp(t)cr z /2)[ax + cr z ] with tan^(t) = 
^Jpy + [A/v F ) 2 l\Vx — eA(t)]. The resulting instanta- 
neous energy spectrum in the upper Dirac cone is e p (t) = 

The transformed time 



A 2 



v 2 F (( Px - eA(t)) 



Pi)- 



dependent Dirac equation reads as 



ihd t $ p (t) 



a z e p (t) 



A 2 



2e 2 Jt) 



(2) 
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FIG. 1: (Color online) The short time Hall conductivity is 
shown from Eq. @ (red dashed line) together with the nu- 
merical solution of Eq. ([2| (blue solid line) for weak electric 
field. The inset shows the numerical results for strong fields 
with the characteristic oscillations set by the field. The tran- 
sient response in both cases is well described by Eq. ([5]). 



and ^p(t) = U$>p(t), with initial (ground-state) condi- 
tion &p(t = 0) = (0,1), in which the lower (upper) 
Dirac cone is fully occupied (empty). The electric held al- 
ters the energy spectrum and induces off-diagonal terms 
in the Hamiltonian. Two energy scales at the moving 
Dirac point (p = (eEt,0)) in Eq. @, characterise the 
low-energy physics: the diagonal energy (A) and off- 
diagonal coupling (hv F eE/2A), which triggers transi- 
tions between the two gap edges or levels using Landau- 
Zener terminology (l7| . A crossover from weak to strong 
field is thus expected at E ~ A 2 /hv F e, irrespective of 
the explicit value of t, as we confirm below by a more 
detailed analysis. 

The quantity we focus on is the time dependent trans- 
verse charge current, j y = —ev F a y in the basis of Eq. 
([T]), with spin current and conductivity differing only by 
a factor h/ev F . For TI, j y coincides with the topologi- 
cal magnetic field induced parallel to the applied electric 
field (after the 7r/2 rotation of the spin, leading to Eq. 
(H])) and monitors the magnetoelectric effect @, 

By denoting $J(f) = (a p (t) , (3 p (t)) , charge conserva- 
tion implies |/3 p (t)| 2 = 1 — n(t), where n p (t) = |a p (t)| 2 
is the number of electrons which have tunneled into the 
initially empty upper Dirac cone. After multiplying the 
transformed Dirac equation, Eq. @ with $+(i)cr x or 
$p(t) from the left, we get 



(jy)p(t) 



ev F AH ( d t [e 2 p (t)d t n p (t)] 
- p {t) yv F eE((v F p y ) 2 

v F eE . , . 

(2n p (f) - 1 



+ A 2 



■+ 



(3) 



which depends only on n p (t) and its time derivatives. 

We start (Fig. [T]) by analyzing its behaviour at weak 
electric fields, determined by the first term in Eq. ([3]). 



The time dependent Hall current is 

h 



^ {t) = 2h 



2 sin" i£ 
h 



E. 
(4) 



where Si(a;) is the sine integral, and exhibits damped 
oscillations around the quantized value of the Hall con- 
ductivity as <T xy = e 2 /2h[l + hsm(2At/h)/2At] with a 
frequency of 2A/h. Still at short times, but in the oppo- 
site small A and strong E limit, similar oscillations with 
a frequency ~ \J v F eE/h show up in the response around 
the non-quantized asymptotic value. The transient be- 
haviour at very short times rises linearly with AEt as 



3y(t) 



irAt 



mm(l, Wt/h)E 



(5) 
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with W the cut-off, and is universal (no high energy 

modes involved) unless t < h/W. 

In the long time limit (t »min[7i/A, ^Jh/v F eE]), we 
can use the analogy of Eq. ([2]) to the Landau-Zcncr 
problem [9I [l7j of two-level crossing to determine n p (t): 



n p {t) = Q(p x (eEt-p x ))eyLp 



TT[(v F p y ) 2 + A 2 

vpheE 



(6) 



which is the pair production rate by Schwinger[ll| and 
also the Landau Zener transition probability |17| between 
the initial and final levels, applicable if ip x , eEt — p x ) 3> 

^JpI + (A/ VF y. 

In the limit of long times, the second term in Eq. ^ 
dominates, and the transverse current reaches a time in- 
dependent value j y (t) = <J xy E with 

{ev F ) 2 A 



>xy 
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which is our main result, eri(x) being the error function. 

Interestingly, the structure of the non-equilibrium Hall 
conductivity at long times agrees with the conventional 
equilibrium Kubo expression [l8[ Ell after shifting the 
momentum with the vector potential and replacing the 
equilibrium Fermi functions with the non-equilibrium 
momentum distribution, Eq. (|6]). Alternatively, it re- 
flects the competition between Berry's curvature (Q p = 
v F A/2e p (t)), protecting quantization and the difference 
of momentum distributions in the upper (n p (t)) and 
lower (1 — n p (t)) Dirac cones, spoiling it: when the two 
distributions are comparable due to tunneling, the gap 
becomes irrelevant, and the conductivity decays. In the 
limit of small fields (E <C irA 2 /vphe), we recover the 
quantized value 



>xy 



e- I Q p e 2 

- ' dp 2^ = Yw 



(8) 



without higher order perturbative or power law (in E) 
corrections. The additional terms contain the non- 
perturbative, exponential factor exp(— irA 2 /vpheE). sig- 
naling the robustness of Hall quantization [20 ] and the 
half-integer quantized magnetoelectric polarizability [7] . 
In the strong field limit (E ^> nA 2 /vF?ie), it decays as 



2A 



2h y/v F heE' 



(9) 



For TI with a mass gap (A ^ 0), the magnetization pro- 
duced by surface currents probes the Hall conductivity 
through the topological magnetoelectric effect, and the 
magnetization parallel to the electric field follows Eq. 
(J7j): its quantization breaks down with increasing field 
similarly to the Hall response. When A = 0, the magne- 
tization perpendicular to E becomes finite ~ (e 2 Tr/2h)E 
in weak fields 

Assuming a small gap of the order of 0.01-1 K (typ- 
ical for the intrinsic SOC of graphene [H [2l| or TI) the 
crossover field is 0.001-10 V/m for vf ~ 10 6 m/s, eas- 
ily accessible experimentally. The Hall conductivity to- 
gether with numerical results on the Dirac equation is 
shown in Fig. [2] The agreement between the analytically 
and numerically obtained conductivities is excellent. 

We can get acquainted with the above results in dif- 
ferent ways: first, a similar situation occurs within equi- 
librium linear response (small E): the (spin-) Hall con- 
ductivity of massive Dirac electrons (i.e. graphene with 
intrinsic SOC and Tl[22j|) is quantized to e 2 /2h, when 
the chemical potential lies within the gap (\fj,\ < A). 
For | fi | > A, the chemical potential cuts into the con- 
tinuum of band states, and the conductivity decays as 
e 2 A/2h\(i\ 1 surviving even the effect of disorder [23[ \24\. 
Within our time-dependent formalism, the electric field 
can be thought formally as introducing an effective chem- 
ical potential /i e // ~ y/HvFeE (only |uf£>j/| < (J-eff con- 
tributes), which upon substituting into the above linear 
response expressions, parallels our findings. 
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FIG. 2: (Color online) The long time limit of the Hall con- 
ductivity is plotted as a function of the applied longitudinal 
electric field. Quantization breaks down when E ~ A 2 /hvFe. 
The red circles denote the numerical data from brute force 
integration of the Dirac equation, Eq. ([2}, while the black 
dashed line is the approximate expression from Eq. ((9} at 
large fields. 



Second, in the edge state picture, the Hall conductiv- 
ity is provided by gapless, one-dimensional ballistic edge 
states, giving rise to the quantized value, which holds for 
weak electric fields. For strong fields, another type of 
gapless excitation starts to contribute, due to tunneling 
between valence and conduction bands (Schwingcr's pair 
production or Zcncr's dielectric breakdown), spoiling the 
perfect Hall quantization, as demonstrated above. 

Another way to look at it is to consider the comple- 
mentary stationary problem to Eq. ([T]) of a static electric 
field in the form of a scalar potential (~ eEx), and an- 
alyze the evolution of the spectrum and edge states as a 
function of the electric field. As a tight-binding example, 
we consider the spectrum of a zigzag graphene ribbonQ 
with intrinsic SOC, causing a gap with opposite sign be- 
tween the two valleys, sublattices and spin directions in 
the continuum limit. In the absence of an electric field, 
only the edge states, connecting the two Dirac cones, 
carry the transverse current, while in a strong electric 
field perpendicular to the edges, the effect of edge states 
is supplemented by the appearance of additional low en- 
ergy modes living in two dimensions, due to the bands 
approaching each other, as seen in Fig. |3] 

At the same time, the longitudinal conductivity is con- 
siderably enhanced in a strong electric field. In the time 
dependent framework, the increasing number of electron- 
hole pairs, tunneled through the gap, facilitate longitu- 
dinal transport while in the stationary picture of Fig. 
I31 the presence of new conducting channels around zero 
energy contribute the electric response. Note that the 
time-dependent framework provides a finite (spin-) Hall 
conductivity even in the absence of scattering [19j,[23j due 
to its intrinsic character, not unlike the analysis of metal- 
lic graphene [2~ij . where additional disorder induced cor- 
rections were found, which we also expect to occur when 
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scattering is added to our framework; these will also limit 
the longitudinal conductivity Q. 




efficient decay as E 1 / 2 . 
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FIG. 3: (Color online) The energy spectrum of a spin-Hall 
insulator[3j] in graphene, t cc is the hopping. Left panel: 
without electric field, showing two gapless, spin degener- 
ate edge states. Right panel: with finite critical electric 
field (red/black denoting up/down spin states), distorting the 
spectrum, and bringing additional levels into play around zero 
energy. Consequently, the spin-Hall conductivity is not quan- 
tized any more. Similar effects are generated by a strain in- 
duced pseudo-electric field having opposite sign in the two 
valleys, resulting in a valley-Hall effect. For stronger E, band 
crossing is more significant. 

Third, 2D Dirac electrons in crossed stationary in 
plane electric (E) and perpendicular magnetic (B) field 
exhibit Landau quantization and subsequently quantized 
Hall conductivity. However, at E = v F B, all Landau 



levels collapse [25l427| and should give way to a different 
Hall response. By defining the energy gap as the distance 
between the Landau levels closest to the Dirac point, we 
get A Landau = v F V2heB, yielding E = A 2 Landau /2Hev F 
for the field causing the collapse of Landau levels, which 
agrees well with the crossover field where the (spin-) Hall 
response changes dramatically. We expect that some of 
our results can be transcribed to the quantum Hall break- 
down in graphene [l^|, testified by a Hall conductivity 
decreasing with the electric field, similarly to Eq. ([9]) . 

These results are relatively robust against disorder, be- 
cause the basic ingredient of the calculation is the non- 
equilibrium momentum distribution function, Eq. ([6]). 
which follows also from a semi-classical approach (WKB, 
expansion in h). 

We have studied the breakdown of the spin-Hall effect 
in graphene and surface Hall and magnetoelectric effect 
in topological insulators at nonzero electric field by con- 
centrating on the real time dynamics of the transverse 
current. The quantization of a xy , protected by the topo- 
logically invariant Chern number [4j, remains intact as 
long as the Hamiltonian varies smoothly (i.e. weak elec- 
tric fields). When non-adiabaticity enters via Landau- 
Zener transitions in strong field, quantization is lost and 
the Hall conductivity as well as the magnetoelectric co- 
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